Draft version May 8, 2012 

Preprint typeset using L£-T^X style cmulatcapj v. 03/07/07 



(N 

o 



X 

Oh' 

(N 

in 

m 
in 



X 



THE PECULIAR EVOLUTIONARY HISTORY OF ICR J17480-2446 IN TERZAN 5 
Alessandro Patruno l , M. Ali Alpar 1,2 , Michiel van der Klis l , Ed P.J. van den Heuvel 1 

Draft version May 8, 2012 

ABSTRACT 

The low mass X-ray binary (LMXB) IGR J17480-2446 in the globular cluster Terzan 5 harbors an 
11 Hz accreting pulsar. This is the first object discovered in a globular cluster with a pulsar spinning 
at such low rate. The accreting pulsar is anomalous because its characteristics are very different from 
the other five known slow accreting pulsars in galactic LMXBs. Many features of the 11 Hz pulsar 
are instead very similar to those of accreting millisecond pulsars, spinning at frequencies > 100Hz. 
Understanding this anomaly is valuable because IGR J17480-2446 can be the only accreting pulsar 
discovered so far which is in the process of becoming an accreting millisecond pulsar. We first verify 
that the neutron star (NS) in IGR J17480-2446 is indeed spinning up by carefully analyzing X-ray 
data with coherent timing techniques that account for the presence of timing noise. We then study 
the present Roche lobe overflow epoch and the two previous spin-down epochs dominated by magneto 
dipole radiation and stellar wind accretion. We find that IGR J17480-2446 is very likely a mildly 
recycled pulsar and suggest that it has started a spin-up phase in an exceptionally recent time, that 
has lasted less than a few 10 7 yr. We also find that the total age of the binary is surprisingly low 
( & 10 8 yr) when considering typical parameters for the newborn NS and propose different scenarios 
to explain this anomaly. 

Subject headings: stars: neutron — X-rays: stars 



1. INTRODUCTION 

According to the recycling scenario (Alpar et al. 
1982, Radhakrishnan & Srinivasan 1982), millisecond ra- 
dio pulsars are produced by spin-up of the neutron stars 
(NSs) via transfer of angular momentum through accre- 
tion in low mass X-ray binaries (LMXBs). During this 
phase, channeled accretion of plasma along the NS mag- 
netic field lines can produce X-ray pulsations that reveal 
the spin frequency of the accreting NS. A confirmation of 
this scenario came with the discovery of the first accret- 
ing millisecond X-ray pulsar SAX J1808. 4-3658, spinning 
at a frequency of about 401 Hz (Wijnands & van der Klis 
1998). Thirteen more accreting millisecond pulsars have 
been discovered so far with spin frequencies ranging from 
about 180 up to 600 Hz. 

In 2010 October the new accreting pulsar IGR J17480- 
2446 has been discovered in the globular cluster Terzan 
5 (Bordas et al. 2010; Markwardt & Strohmayer 2010). 
Pooley et al. (2010) localized the source in outburst with 
Chandra observations and identified it with a previously 
known candidate quiescent LMXB, named CX25 in the 
X-ray source catalog of Heinke et al. (2006). The source 
is in core of Terzan 5, one of the densest among globular 
clusters (Ferraro et al. 2009). 

The accreting pulsar spins with a frequency of 11 Hz 
and the LMXB has a period of 21.3 hr around a com- 
panion with mass M > 0.4 M Q (Papitto et al. 2011). 
Although this pulsar is not a millisecond one, its loca- 
tion in a globular cluster makes it a unique system whose 
properties are of particular importance to study the re- 
cycling mechanism. IGR J 17480-2446 is different from 
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the other five known slow accreting pulsars in LMXBs 
(4U 1626-67, 2A 1822-371, GRO 1744-28, GX 1+4 and 
Her X-l) with properties which are in between those sys- 
tems and the accreting millisecond pulsars (Linares et al. 
2012). Slow accreting pulsars are a very heterogeneous 
group of LMXBs, with white dwarf, main sequence and 
giant companions and orbital periods ranging from less 
than 1 hr to more than 11 days. Despite this, the five 
slow pulsars have a high magnetic field B ~ 10 11 — 10 13 G 
and do not show thermonuclear bursts. This high B field 
is difficult to reconcile with the long lifetime of these slow 
systems, and several scenarios have been proposed (see, 
for example, Verbunt et al. 1990 for a discussion). 

IGR J17480-2446 instead shows thermonuclear bursts 
(Chenevez et al. 2010, Linares et al. 2011b) with burst 
oscillations phase locked with the accretion powered pul- 
sations (Cavecchi et al. 2011), and the B field has been 
constrained to be in the range 10 9 — 10 10 G (Miller et al. 
2011, Cavecchi et al. 2011, Papitto et al. 2011). All 
these field estimates point toward a mildly recycled pul- 
sar, since accreting millisecond pulsars have shown so far 
field strengths of the order of 10 8 G (in agreement with 
the B fields of radio millisecond pulsars). Phase locking 
between burst oscillations and accretion powered pulsa- 
tions has also been observed in at least one accreting 
millisecond pulsar (Strohmayer et al. 2003, Watts et al. 
2008). The 11 Hz rotation rate is, however , still unusu- 
ally slow for an old NS, the typical rotation periods for 
accreting or rotation powered pulsars in globular clus- 
ters and old stellar populations being in the millisecond 
range. 

The reason for old NSs possessing dipole magnetic 
fields smaller by about three orders of magnitude from 
the typical 10 12 G fields of young NSs has elicited several 
explanations involving binary evolution. The burial of 
the field under accreted material has been proposed as a 
viable mechanism to dissipate the field via ohmic decay 
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in the heated NS crust (see for example Romani 1990). 
This mechanism, however, has a dramatic consequence 
for the final magnetic field: as the accreted matter keeps 
falling on the NS surface, the magnetic field lines will 
be pushed to deeper (and colder) regions of the crust 
where the conductivity is much larger, thus freezing at 
the crust-core interface. This will leave a residual field 
that is compatible with the observed weak fields of mil- 
lisecond pulsars. 

Another particularly simple and elegant explanation 
relies on the superfluid-superconducting nature of the NS 
interior, requiring that the rotation is carried by quan- 
tized vortex lines in the neutron superfluid and the mag- 
netic flux by quantized flux lines in the proton super- 
conductor (Srinivasan et al. 1990). Spin down of the 
star is achieved by outward motion of the vortex lines, 
which inevitably entangle and drag flux lines with them. 
Srinivasan et al. noted that many NS binary evolution 
scenarios have an early first mass transfer epoch when 
the companion, not yet filling its Roche lobe, is losing 
mass by a stellar wind, some of which is captured by the 
NS in a weak quasi spherical accretion with low specific 
angular momentum. The NS will be spinning down in 
this epoch, thus reducing its dipole field in proportion to 
its rotation rate. The wind spin-down epoch ends when 
the companion fills its Roche lobe. The ensuing disk ac- 
cretion now starts to spin- up the NS. 

If IGR J 17480-2446 is a primordial binary, then its spin 
frequency is surprising low (the globular cluster Terzan 5 
is about 12 Gyr old, see Ferraro et al. 2009) and requires 
investigation to understand whether this is the conse- 
quence of a peculiar evolutionary history. At present IGR 
J 17480-2446 is likely to be in its spin- up epoch (Papitto 
et al. 2011). The results reported by Papitto et al. (2011) 
are however affected by the usual problem of X-ray tim- 
ing noise in the pulse time of arrivals (TOAs) that when 
not taken properly into account might strongly affect the 
significance of the spin frequency derivatives (see for ex- 
ample the case of the other three accreting millisecond 
X-ray pulsars in Hartman et al. 2008; Hartman et al. 
2009, Patruno et al. 2010 Haskell & Patruno 2011). In 
this work we will re- investigate the existence of a spin-up 
phase by taking into account the effect of timing noise. 
If the spin-up is confirmed, then the comparatively long 
spin period indicates either that IGR J17480-2446 is at 
an exceptionally early phase of the spin-up epoch or that 
it is close to spin equilibrium, with an unusually long 
equilibrium period reflecting a magnetic field stronger 
than in most NSs in LMXB. 

The aperiodic variability can also help to under- 
stand the peculiar behavior of the source. Altamirano 
et al. (2010) reported several quasi-periodic oscillations 
(QPOs) at frequencies ranging from about 48 up to 815 
Hz. This sequence of frequencies make possible sev- 
eral alternative estimations of the inner disk radius at 
present. Each such model gives an estimate of the dipole 
magnetic field, equilibrium rotation rate and expected 
spin- up rate. 

The location of IGR J17480-2446 in the globular clus- 
ter Terzan 5 allows also to place stronger constraints on 
the evolutionary history of the binary since only a nar- 
row range of donor masses is allowed. The distance of 
the system is also constrained by the globular cluster and 
leads to a precise determination of the X-ray luminosity. 



In Section 2 we describe the X-ray observations carried 
during the 2010 outburst, and in Section 3 we verify the 
existence of a spin-up phase. In Section 4 we present the 
evolutionary history of IGR J 17480-2446 . We discuss 
three distinct epochs: in the initial two the NS has been 
spun down first via magneto dipole radiation and then 
by wind accretion, whereas in the third epoch the NS is 
spun-up via accretion of mass provided by the donor in 
Roche lobe overflow (RLOF). The evolution of the mag- 
netic moment of the NS through the earlier dipole and 
wind spin-down epochs is also discussed. In Section 5 
we place further constraints on the binary parameters 
by using the available information on the globular clus- 
ter Terzan 5 and the results of the evolutionary history. 
We conclude with a discussion of the results in Section 6 
and infer that IGR J 17480-2446 is in an exceptionally 
early RLOF phase that explains why its spin frequency 
is still so small when compared to accreting millisecond 
pulsars. 

2. X-RAY OBSERVATIONS 

We reduced all pointed X-ray observations taken with 
the Rossi X-ray Timing Explorer (RXTE) by the Pro- 
portional Camera Array (Jahoda et al. 2006) in event 
mode. Only photons with energies between 2.5 and 15.6 
keV are retained (absolute channels 5—37). The time 
resolution of the observations is 2~ 13 s (for 122/j,s Event 
data) and 2 -20 s (for Good Xenon data) rebinned to 2~ 13 
s to make the data sampling uniform. The photons de- 
tected span an observing window between 2010 October 
13 and November 19. The total duration of the outburst 
is longer than the duration of the RXTE observations, 
and it has been constrained to be « 55 days (Cavecchi 
et al. 2011) or even longer (77 days) when modeling the 
outburst decay (Degenaar & Wijnands 2011). 

2.1. Data Reduction Procedure for Coherent Analysis 

A previous coherent timing analysis of IGR J17480- 
2446 is briefly described in Cavecchi et al. (2011), 
who reported a pulse frequency derivative of « 1.4 x 
10 _12 Hzs _1 with errors calculated by means of Monte 
Carlo (MC) simulations that take into account the red 
noise in the timing residuals (Hartman et al. 2008). The 
steps followed are those used in standard coherent timing 
of radio and X-ray pulsars: the photon arrival times are 
referred to the solar system barycenter using the optical 
position from Pooley et al. (2010), fine clock corrections 
are applied to the RXTE data, and Earth occultations 
and intervals of unstable pointing are filtered out. 

To measure the TOAs of the persistent pulsations, all 
bursts are removed from the light-curve, then w 10 to 
« 500 s long intervals of data are folded in profiles of 
N = 32 bins with the preliminary ephemeris reported in 
Papitto et al. (2011). The length of the folded data seg- 
ments is chosen to have a nearly uniform statistical error 
(~ 1ms) for each TOA. The « 10 s profiles are those 
folded from the first RXTE data-set (Obsld 95437-01- 
01-00) when the pulse profiles have the highest fractional 
amplitude and the highest signal-to-noise ratio. All other 
observations are folded with w 500 s segments with the 
exact length slightly varying for each observation to as- 
sure that all good photons are retained in the data anal- 
ysis. We assume that pulsations are detected in a profile 
if the ratio between the amplitude of a harmonic and 
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its statistical error is larger than 3.5. This guarantees 
that the expected number of false pulse detections is less 
than one in the entire data collection (the total number 
of profiles is w 650 for each harmonic). A fundamental 
and two overtones are detected with strong significance 
and in enough data segments to allow a determination of 
a precise global timing solution. 

The TOAs of each harmonic are fitted using the 
TEMP 02 pulse timing program (Hobbs et al. 2006), as- 
suming a circular Keplerian orbit plus a pulse frequency 
and its first time derivative. The entire folding and fit- 
ting procedure is then repeated until the final timing so- 
lution converges. The orbital parameters are reported in 
Table 1 along with 1 a error bars obtained by multiply- 
ing the statistical error as given by TEMP02 by \fx*- 
The errors on the pulse frequency and its derivative are 
instead treated separately to include the strong effect 
of timing noise which operates on the same timescales 
over which the two parameters are measured. Details 
of the procedure have been extensively discussed in the 
literature, and we refer to Cavecchi et al. (2011) for the 
results obtained on IGR J17480-2446 and to Hartman 
et al. (2008) and Patruno et al. (2010) for further details 
on the procedure. 

Here we will verify more carefully whether the pulse 
frequency derivative (i.e., the measured i>) reported in 
Cavecchi et al. (2011) is indeed the spin frequency deriva- 
tive (i.e., a true rotational variation of the NS spin) and 
not the effect of timing noise. Although Cavecchi et al. 
(2011) used MC simulations to calculate the long-term 
average v of IGR J 17480-2446 , no interpretation to the 
pulse frequency derivative is given and no discussion of 
the effect of timing noise is provided. Another long-term 
v has been reported by Papitto et al. (2011). It is of the 
same order as, but significantly different from, that given 
in Cavecchi et al. (2011). The timing analysis of Papitto 
et al. (2011) did not consider in detail the strong effect of 
timing noise and the authors reported discordant v val- 
ues in their analysis of the fundamental and first overtone 
which deviate by more than 32ct from each other and de- 
viate by more than 3cr from the results of Cavecchi et al. 
(2011). Identical considerations also apply to the pulse 
frequency. These results show why extra care has to be 
taken to verify that the pulse frequency derivative is the 
spin frequency derivative of IGR J 17480-2446 . Deter- 
mining whether the NS is in a spin-up phase is crucial 
for understanding the evolutionary history of the binary 
as we will discuss in Section 4. 

We split the data in five segments of length between 
« 4 and 9 days (see Table 2). The exact values are cho- 
sen empirically to have enough data to measure a pulse 
frequency derivative of the order of 10~ 12 Hzs _1 . The 
data used and the start and end times of each data seg- 
ment are reported in Table 2. We fit the pulse frequency 
and its first derivative to each data segment. The fun- 
damental frequency and the first overtone are analyzed 
separately, whereas the second overtone is not detected 
in enough observations to allow data splitting and will 
therefore not be consider any further. The harmonic 
phase residuals relative to our best constant frequency 
timing model exhibit variability on time scales as short as 
a few tens of minutes (i.e., close to the shortest timescale 
we are sensitive to) with an amplitude well in excess of 



the Poisson noise expected from counting statistics. This 
variability requires a special treatment of the statistical 
errors of the pulsar spin parameters since standard x 2 
minimization techniques (i.e., those used by TEMP02 or 
any other standard fitting routine) rely on the assump- 
tion that only Poissonian noise is present in the data, 
which is not a good approximation in IGR J17480-2446 
. We apply the MC method developed by Hartman et al. 
(2008) to estimate the effect of timing noise on the pulse 
parameters. We run 10 4 MC simulations on time se- 
ries which have the same power content of the original 
TOA residuals. The parameters v and v are measured 
for each simulated time series and a distribution of pa- 
rameters is constructed over the 10 4 simulations. The 
68% confidence interval on the pulse parameters is then 
determined as the standard deviation of each parameter 
distribution. 

A second test is performed by repeating the entire pro- 
cedure described above with the difference that we fit 
only a constant pulse frequency with the derivative set 
to zero. In this way we perform a "partial coherent" 
analysis: coherent timing is used in each of the five data 
segments to extract a constant pulse frequency (constant 
in that segment), whereas the long term pulse frequency 
derivative is inferred later by fitting the five pulse fre- 
quencies with standard x 2 minimization techniques. In 
this way we can check whether the five constant pulse 
frequencies increase linearly in time as expected if a con- 
stant spin frequency derivative is present in the data. If 
instead the five frequencies fluctuate with a large scatter 
around a linear trend, the effect of timing noise is pre- 
dominant and a presence of a constant spin frequency 
derivative is questionable. 

To establish whether the pulse frequency derivative is 
consistent with a spin frequency derivative we require 
that: 

• the pulse frequency derivative inferred for the two 
harmonics is consistent with being the same and it 
also matches with the long-term pulse frequency 
derivatives measured with coherent timing (i.e., 
those published in Table 1 by Cavecchi et al. 2011); 

• the pulse frequencies have to be consistent with 
each other for the two harmonics in each of the five 
segments; 

• there is no large deviation (> 3tr) of any of the 
pulse frequencies from the linear trend. This 
demonstrates that timing noise is not dominating 
the fitted pulse frequency parameter in any of the 
five segments. 

3. RESULTS OF COHERENT ANALYSIS 

We produce two different outputs: the pulse frequency 
derivative as it evolves in time (Figure 1, left panel) and 
as a function of the X-ray flux (Figure 1, right panel). 
The parameter shows no clear evolution with time and 
the distribution of points around the long-term average 
v is consistent with random fluctuations due to mea- 
surement errors, with the exception of one outlier in the 
fourth data segment of the fundamental. 

Under the assumption that the X-ray flux tracks the 
mass accretion rate M and that our five measured pulse 
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TABLE 1 

Orbital Solution of IGR J 17480-2446 



Orbital Parameter 


Fundamental 


First Overtone 


Second Overtone 


Pb [hr] 


21.27458(7) 


21.27470(7) 


21.2746(2) 


a x sin(i) [lt-s] 


2.4965(4) 


2.4962(4) 


2.498(2) 


T asc [MJD] 


55481.78025(4) 


55481.78015(5) 


55481.78018(2) 


e [95% c.l.] 


< 0.002 


< 0.004 


< 0.004 



TABLE 2 
X-Ray Observations 



Segment Nr. 


Obs Ids (95437) 


Start 
(MJD) 


End 
(MJD) 


Flux 

(Crab) 


1 


01-01-* 01-02-* 
01-03-* 01-04-* 01-05-* 


55482.006 


55486.755 


0.411 


2 


01-06-* 01-07-* 01-08-* 
01-09-* 01-10-00 01-10-01 
01-10-02 


55487.306 


55493.347 


0.499 


3 


01-10-03 01-10-04 01-10-05 
01-10-06 01-10-07 01-11-03 


55494.246 


55501.326 


0.382 


4 


01-11-04 01-11-05 01-11-06 
01-11-08 01-12-* 


55502.407 


55511.115 


0.316 


5 


01-13-* 01-14-* 


55512.153 


55519.165 


0.259 
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Fig. 1. — Left Panel: Time evolution of the pulse frequency derivative of IGR J17480-2446 for the fundamental (blue squares) and first 
overtone (red circles). The dotted line refers to the long-term average v of the fundamental reported in Cavecchi et al. (2011), with the 
green band being the 68% confidence interval calculated with MC simulations. Each data point is measured in a ~ 4 — 9d data interval and 
the errors are calculated by means of MC simulations (see Section 2). All the data points are consistent with the green confidence interval 
bar with the exception of the 4th data point of the fundamental which can be considered an outlier (see Section 3.1 for a discussion). Right 
Panel: Pulse frequency derivative vs. X-ray flux. The meaning of symbols is the same as in the left panel. No trend is visible in the data 
suggesting that short-term v measurements have a very weak or no dependence on the X-ray flux. 



frequency derivatives represent short-term instantaneous 
spin frequency derivatives, we expect to see i> oc f x , 
where 7 = | in the simplest accretion model (see for ex- 
ample Ghosh & Lamb 1979, Bildsten et al. 1997). From 
the right panel of Figure 1 we see no clear dependence of 
v from fx ■ All short term pulse frequency derivatives are 
consistent within 1 a with the average long-term pulse 
frequency derivative v = 1.42(5) x 10 -12 , reported by 
Cavecchi et al. (2011), except the second data point of 
the fundamental frequency which deviates by more than 
3(7 and is the same outlier seen in the left panel of Fig- 
ure 1 (see Section 3.1 for further discussion of the outlier). 
Finally, in Figure 2 the constant pulse frequencies are 



plotted against time for fundamental and first overtone. 
We fit each harmonic with a linear function f(t) = k + i>t 
to extract the long-term pulse frequency derivative. The 
results give i> and k for the two harmonics which are 
consistent within the la statistical errors and are both 
consistent within la with the long-term pulse frequency 
derivative reported in Cavecchi et al. (2011). 

3.1. The Spin-up of IGR J17480-2446 

The results presented above strongly suggest that the 
long-term pulse frequency derivative reported by Cavec- 
chi et al. (2011) is indeed the long term spin frequency 
derivative of IGR J17480-2446 . There are however two 
points that need further discussion before drawing a firm 
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Fig. 2. — Pulse frequency evolution with time. The data points 
refer to the same data segments used in Figure 1 and described in 
the text. The blue squares and red circles refer to fundamental and 
first overtone. Each pulse frequency is calculated with a constant 
pulse frequency model. Both data-sets (fundamental and first over- 
tone) can be fit with a linear relation (blue dotted: fundamental, 
red dashed: first overtone) which is consistent with being the same 
within the statistical errors and gives a long-term v consistent with 
the value reported in Cavecchi et al. (2011). 

conclusion from our data analysis: the meaning of the 
outlier in the short-term pulse frequency derivatives and 
the lack of correlation between v and fx ■ 

Although the outlier indicates a significantly larger 
pulse frequency derivative than the long-term average, 
the first overtone gives in the same data segment a v 
incompatible with the outlier and compatible with the 
long-term average v. The interpretation we give to this 
is that the outlier does not represent a true increase in 
the spin frequency derivative of the NS. This can be still 
ascribed to timing noise, even if we are carefully calculat- 
ing the statistical errors with MC simulations. MC sim- 
ulations verify whether a measured v is significant with 
respect to the red noise component in the TOA residuals. 
This is achieved by testing whether random fluctuations 
in the residuals are able to produce fake v as strong as 
the measured one. Since we are looking at short time 
series (~ 4 — 9days), the red noise is not stationary, as 
evident from Figure 3, where the data segment to which 
the outlier belongs is shaded. In that data segment, the 
highest red noise frequencies observable are comparable 
with the length of the data segment and the data show a 
simple parabolic trend which is absorbed by v during the 
fitting procedure. If this is the case, MC simulations can- 
not discriminate between a true spin frequency derivative 
and the effect of timing noise. Therefore the error bars 
thus calculated do not differ much from those one obtains 
when the errors are dominated by uncorrelated Gaussian 
noise (i.e., white noise). The final values of the MC er- 
rors are not, is such good representation of the 
true statistical uncertainties. This is a limitation of the 
MC method. However, a similar behavior is not observed 
in the same data segment for the first overtone and so 
we ascribe this behavior to timing noise. 

It is of course possible that timing noise is the result of 
true changes of the spin frequency and / or spin frequency 
derivative of the NS. Fluctuating torques have been pro- 
posed in the past as a mechanism to explain fluctuations 
in the NS rotation (Lamb et al. 1978). However, in this 
case we would expect to see identical variations of the 



spin parameters in both harmonics which is not observed 
for the outlier in ICR J 17480-2446 . Furthermore it has 
been shown that the timing noise observed in at least 
two accreting millisecond pulsars (Watts et al. 2008; Pa- 
truno et al. 2010) cannot be ascribed to true variations 
of the neutron star spin parameters since these require 
unphysical large torque variations. 

The fact that v is not proportional to fx can be a con- 
sequence of different phenomena. The X-ray flux might 
be not a good tracer of M as has been already seen in 
several other X-ray binaries (van der Klis 2001). Our 
data suggest that v k. const over an X-ray flux variation 
of a factor w 2. An approximately constant bolometric 
flux might therefore easily explain the observations. 

Another possibility is that v oc /j^ 7 is valid only if the 
accretion disk is a standard thin disk under the rather 
unrealistic assumption that radiation pressure is not im- 
portant in determining the inner disk structure (Psaltis 
& Chakrabarty 1999). Andersson et al. (2005) suggest 
that this assumption is not correct already at accretion 
rates of a few percent Eddington. Since IGR J17480- 
2446 is accreting at least at M > 5% Eddington with 
peaks exceeding 37%, it is not unreasonable to assume 
that radiation pressure plays a role in determining the in- 
ner disk structure. The disk-magnctosphere interaction 
region is also more complex than usually assumed (see 
e.g., Rappaport et al. 2004, D'Angelo & Spruit 2010, Ka- 
java et al. 2011, Patruno et al. 2012), and v might have 
a weaker dependence on the bolometric flux than usually 
assumed. 

We conclude therefore that the interpretation of the 
long term v as being the spin frequency derivative of 
IGR J17480-2446 is supported by the observations, with 
the outlier most likely explained by unmodeled timing 
noise and the lack of correlation between the short term 
v and fx, which is not surprising, since the expectation 
of such correlations is based on assumptions about the 
accretion process which may not be valid. 

4. EVOLUTIONARY HISTORY OF IGR J17480-2446 

The present RLOF spin-up epoch has been preceded 
by two earlier epochs in which the binary was still de- 
tached. In the first dipole spin-down epoch, the newly 
born NS has an initial rotation rate fio (spin vq = flo/2ir) 
and dipole moment fi ■ The magneto dipole torque dom- 
inates this epoch until the wind of the companion pen- 
etrates the light cylinder and the NS is spun down by 
wind torques. When this second epoch starts, the NS 
has rotational rate and dipole moments 51i and fii re- 
spectively. In this section we present an evolutionary 
history of these three epochs to determine the evolution 
of the NS spin in IGR J17480-2446 . The dipole mo- 
ment decay induced by the spin-down of the NS in the 
superfluid-superconducting core of the star is assumed to 
be an essential factor during the spin-down epochs. We 
summarize the duration and the symbols used for each 
evolutionary phase in Table 3. 

4.1. Duration of the Roche Lobe Overflow Epoch 

In the present RLOF spin-up epoch, the 2-16 keV 
X-ray luminosity of IGR J17480-2446 spans values be- 
tween 9 x 10 36 ergs~ 1 and 6.5 x 10 37 ergs _1 for an as- 
sumed distance of 5.5 kpc (Ortolani et al. 2007). This 
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Fig. 3. — Timing residuals for the fundamental (left panel) and first overtone (right panel). The strength of timing noise for the two 
harmonics is comparable with an rms of the residuals for fundamental and first overtone of 5.1 and 4.9 ms respectively. The fundamental 
shows a stronger timing noise than the first overtone on short timescales (~ Id) for the first 15 days and the last 3 days, whereas the first 
overtone has a stronger noise component on longer timescales. The yellow shaded area refer to the data segment where the fundamental 
pulse frequency derivative shows an outlier (see text for a discussion) 



TABLE 3 

Evolutionary History of IGR J 17480-2446 



Epoch Duration (yr) /jj fir Qj Of 

Dipolc spin-down 10' JIq H-i = Ho f2o 

Wind spin-down 10 7 — 10 8 fj,\ = [iq fi2 f^i ^2 

RLOF spin-up 10 7 - 10 8 fj® M2 ^2 ^3 



determines the torques on the NS applied by the ac- 
cretion flow and the induced distortions of the magnetic 
field. These torques operate as long as ft is not equal 
to njf(r>i), so that the star spins up or down toward an 
equilibrium rotation rate Q eq = CIk{ta)- If any of the 
QPO frequencies (see Altamirano et al. 2010) is inter- 
preted as the eventual spin equilibrium frequency of the 
NS, the present rotation rate will give an estimate of the 
dipole moment of the NS. This identification however, is 
not straightforward. QPO frequencies reflect modes of 
oscillation in the accretion flow that are excited in the 
interacting disk-NS system to some amplitude to effect 
the luminosity signal. A particular QPO frequency band 
might be defined by some oscillation mode of the disk 
or by a beat or resonance between the NS and a disk 
mode. Disk modes are fundamentally effected by rota- 
tion and imprinted by the frequency scale of the disk 
rotation, Keplerian in most parts of the disk, but nec- 
essarily and significantly deviating from the Keplerian 
rate at the inner boundary where the effect of the NS 
magnetosphere is important. This results in frequen- 
cies of epicyclic (radial) disk oscillations to be somewhat 
higher than the non-Keplerian rotation frequency at the 
inner disk boundary or transition region (Alpar & Psaltis 
2008). The presence of a sequence of QPO frequencies 
may be due to the presence of several multipole com- 
ponents of the neutron star magnetic field, so that the 
boundary between the disk and the NS magnetosphere 
is convoluted with a sequence of decreasing "stopping- 
radii" , smaller than the conventional Alfven radius given 
in Equation (2) and increasing Kepler frequencies corre- 
sponding to higher multipoles of the NS magnetic field 
(Alpar 2011). We assume that one of the observed range 
of QPO frequencies z/j is close to the Keplerian rotation 
rate SIk(ta) near the Alfven radius, 

(4) 



Note. — The subscripts i and / refer to the initial and final 
value at a specific epoch. The magnetic moment 112 is assumed to 
remain constant during the RLOF phase. 

luminosity has to be considered a lower limit given that 
the bolometric luminosity is expected to be higher than 
the X-ray luminosity. Linares et al. (2012) estimated 
a bolometric luminosity at the peak of the outburst 
of about 50% the Eddington limit for d = 6.3 kpc 
(L>Edd — 2.5 x 10 38 ergs _1 ). The maximum mass ac- 
cretion rate is therefore 

i/ = _^vs ^ 5xl0 i 7 f d_Y*L -i (1) 
GM NS V5.5kpcy Mia 

where G is the gravitational constant, Mns and Rns 
are the mass and radius of the NS, and Mi. 4 and Rq are 
the mass of the NS in units of 1.4 Mq and the radius in 
units of 10 6 cm respectively. The average mass accretion 
rate during the outburst has a value approximately half 
the peak value, and we assume it is ( M ) — 2 x 10 17 g s~ 1 , 
in agreement with the value reported by Degenaar & 
Wijnands (2011). 

The inner radius of the accretion disk is determined 
by the balance between the magnetic stresses in the NS 
magnetosphere and the material stresses of the accretion 
disk (Ghosh & Lamb 1979): 

ta= (jQ 1 (GM NS y 1/7 . (2) 

where \x is the dipole magnetic moment of the star. 
The Keplerian rotation rate at the inner radius of the 
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TABLE 4 

QPO Frequency Model for IGR J17480-2446 



QPO Frequency 




B Field 


V 


tsu,QPO 


(Hz) 


(10 27 Gcm 3 ) 


(10 9 G) 


(io- 12 ) 


(10 7 yr) 


48 


5 


10 


1.6 


1.8 


173 


1 


2 


1.1 


2.8 


814 


0.2 


0.4 


0.6 


4.6 



Note. — The magnetic field B refers to the value at the NS 
magnetic poles. This is twice the value at the equator. 

yielding the estimates 



H = (27r^)- 7/6 M 1/2 (GM NS ) 5/e 



(5) 



1.3 x 10 26 G cm 



VkHz/ 



for the dipole moment. In these estimates we use (M) to 
mean the mass accretion rate averaged over the outburst 
length and it is normalized to 2 x 10 17 g s _1 ; 

The values of the magnetic dipole moment are reported 
in Table 4. If we take the NS to be already in equilibrium 
at the present spin frequency of 11 Hz, then \i = 2 x 
10 28 G cm 3 and a B field at the poles of 4 x 10 10 G where 
we have used: 



M30 



0.6/ 7 / 6 



d 



5.5 kpc 



V^M^V 7 / 6 (6) 



Here / takes account of the possibility that the observed 
frequency is not the Kepler frequency but some related 
frequency at the disk boundary: e.g., / ~ 1.5 — 2 for the 
epicyclic frequency (e.g., Alpar & Psaltis 2008). This 
uncertainty makes the field estimates somewhat higher. 
The torque applied by the accretion disk on the NS is 
usually assumed to be: 



N = M(GM NS r A ) 1/2 F(uj) 



(7) 



with the dimensionless torque F(u>) describing the de- 
pendence on the fastness parameter u> = fl /ClK(r A ), is 
likely of the form 



F(u) 



(8) 



with the exponent p ~ 2 (Ertan ct al. 2009). If IGR 
J 17480-2446 is as yet at the beginning of its spin- up 
epoch, then to « 1, F(uj) = 1, giving constant v and 
expected spin frequency derivatives 



i)^6x 10~ 13 Hz s 



-1 (M)M 2/ 4 3 

J 45 (infe) 



1/3 



(9) 



where I45 is the NS moment of inertia in units of 
10 45 g cm 2 . These estimates of v are also given in Table 4. 
With any of the QPO Vi this gives v ~ 10~ 12 Hz s _1 . 
If the system is already in equilibrium v would be 
fluctuating in sign with rms values much smaller than 
10~ 12 Hz s _1 , which is however not observed (see Sec- 
tion 3 and Figures 1 and 2). 

The spin-up age up to the present is then 



1/3 



tsu.gpo — 5 x 10 yr 



7 45(fc7fe) 



(10) 



where A is the (unknown) duty cycle of the binary. If we 
assume A ~ 0.01, the maximum estimated time in the 
spin-up epoch up to now is m 5 x 10 7 yr corresponding to 



Vi = 814 Hz (see Table 4). If any of the QPO frequencies 
is interpreted as the eventual equilibrium frequency, then 
the time spent in the current spin-up epoch is a small 
fraction of the age of the system. Even if the system 
were already in equilibrium at v = 11 Hz then it must 
have reached equilibrium in the short time 



t su = — ~ 2.5 x 10V. 
v A 



(11) 



4.2. 



Duration of the Dipole Spin- down Epoch 

New born NSs have a distribution of rotation rates 
r2o ^ 1 — 300 rad s _1 and dipole magnetic moments 
~ 10 29 -10 30 G cm 3 (Faucher-Giguere & Kaspi 2006). 
The NS is active as a pulsar until the magnetospheric 
voltage falls below a critical threshold (the "pulsar death 
line" ). For conventional radio pulsars the epoch of pulsar 
activity lasts for 10 6 -10 7 yr. No evidence of magnetic 
field decay is seen in the pulsar population during this 
active epoch. The lower limit for an exponential decay 
time is 10 7 yr. 

If the magnetic dipole moment decayed in proportion 
to the rotation rate, due to flux-line vortex line coupling 
(Srinivasan et al. 1990), 

n(t) 



fj,(t) = Mo- 



n 



(12) 



the spin-down would follow a power law with braking 
index = 5. More detailed considerations of flux-line mo- 
tion indicates that the simple scaling in Equation (12) is 
a likely to be a good description for the evolution ( Jahan- 
Miri 2000). 

The vortex- flux coupling leads to the decay of magnetic 
flux in the NS superfluid neutron-superconducting pro- 
ton core. For the total dipole moment to decay, the flux 
expelled from the core must diffuse through the highly 
conducting crust. The decay time through the crust is 
estimated to be 10 7 yr. This is supported by the limits 
on exponential decay of the dipole moment from popula- 
tion analysis of the radio pulsars. Ongoing dipole spin- 
down on time scales much longer than 10 7 yr, that is, 
beyond the typical duration of the epoch of radio pulsar 
emission, will be effected by the power law decay of the 
dipole moment due to flux- vortex coupling. 

Whether pulsar activity has ceased or not, vacuum 
dipole spin-down will continue until the wind captured 
from the companion penetrates the light cylinder. The 
captured wind loses angular momentum in a shock and 
will accrete towards the NS down to a stopping radius 
~ ta- The wind will be able to affect the torque on the 
neutron star when 



r A {p.o,M wNS )^r LC = — • 



(13) 



with M w ns denoting the wind captured by the NS. After 
this point, spin-down will proceed under the effect of the 
inflowing wind. From Equation (13) the rotation rate fii 
at the beginning of the wind spin-down epoch is 

(14) 



fii = 24 rad s^^^M 217 



M 1/7 



where M w ns,13 is the mass inflow rate arriving at ta in 
units of 10 13 g s _1 . If we first assume that flux- vortex 
coupling has no effect, the dipole spin-down proceeds 
with a constant dipole moment and braking index 3: 



3c 3 / 



(15) 



8 



' 2.5 x 10 V n \ 



• -6/7^-4/7 



where c is the speed of light and / is the NS moment 
of inertia. The duration of the dipole spin-down epoch, 
until the rotation rate fii is reached is given by 

. _Pi_ 

" =3 "2P" 2fi 

= 10 1Q yr [i- 2 29 n- 2 

:1' 7 (i6) 

where we have used Equations (14) and (15) and I = 
10 45 gcm 2 . In the last line of Equation (16) we have 
assumed that the magnetic field does not decay during 
the dipole spin-down epoch, so that the dipole moment 
fii at the end of this epoch is still the dipole moment 
fiQ of the NS at birth. Since the lifetime thus obtained 
is not much smaller than the crust diffusion timescale 
of 10 7 yr, it might be necessary to take the vortex-flux 
coupling and the induced decay of the dipole magnetic 
moment into account. The correct estimate of dipole 
spin-down epoch, incorporating flux decay according to 
Equation (12) proceeds with braking index 5. We obtain 



~5x 10 9 yr 4 

= 1.5 x10 V A'LKnIis 



M~ 



4/7 



(17) 

Relating the magnetic moment Hi and rotation rate Q\ 
at the end of dipole spin-down to the initial values Ho 
and f2 with Equation (12) 



in Ql 

"0 



Ml, 29 

and using Equation (14) we obtain 

n 



i ^r, 1n 4 2/11 

t n =5 = 2x10 y /x ; 30 



10 rads 



2/11 



M, 



-12/11 

wNS,13 



(18) 



-6/11 



This value is shorter than the typical field diffusion 
time scale through the NS crust ( k, 10 7 yr) and there- 
fore Equation (19) cannot be used self-consistently for 
Ho ~ 10 29 — 10 30 Gcm 3 . If one ignores the field decay 
and assumes a typical Ho = 10 29 ~ 30 G cm 3 in Equation 
(16) then the dipole spin down age is ~ 10 7 yr. 

4.3. Duration of the Wind Accretion Spin-down Epoch 

Wind accretion (Edgar 2004, Theuns & Jorissen 1993; 
Theuns et al. 1996, Pfahl et al. 2002) proceeds by Bondi- 
Hoyle like capture of material from the companion's 
wind, which is then channeled toward the NS by outward 
transport of angular momentum in traversing a shock 
front. The wind mass loss rate from a ~ 1 M Q main se- 
quence companion is M w ind ~ 10 12 gs" 1 , which is low 
when compared to the present day mass accretion rate, 
driven by RLOF. Even considering a sub-giant phase of a 
~ 1 Mq donor star, the wind loss rate will increase only 
slightly, reaching M w i n d ~ 10 13 gs _1 for sub-giant stars 
(see for example Reimers 1975, 1978). Only a fraction of 
this wind will be captured to flow towards the NS, with 
M/M wind as low as ~ 10 -2 (Nagae et al. 2004). How- 
ever, since the orbit of IGR J17480-2446 shows strin- 
gent upper limits on the eccentricity (e < 10~ 3 , Pa- 
pitto et al. 2011) tidal circularization has probably al- 
ready taken place and it is likely that the donor star is 



rotating synchronously with the orbit, since the circu- 
larization timescale is typically larger than the synchro- 
nization timescale (see for example Hurley et al. 2002). 
The large rotation of the donor (v rot — 70kms _1 ) can 
boost the wind loss rate by a large factor of the order of 
10 2 ~ 3 when compared with typical isolated stars of sim- 
ilar mass at the same evolutionary stage. The captured 
wind might be therefore comparably higher, with values 

ofM^s-lO^-^gs" 1 . 

The NS is spinning rapidly within this wind of low spe- 
cific angular momentum. For our purposes of estimating 
the duration of the spin-down epoch it will be enough to 
take a constant torque, giving the spin-down rate: 



tl = -^{GM NS r A f 2 

t 6/7 M 2/7 ^ ,3/7 



= - M, 



iwNS 



(GM NS y 

-l,V6/7 ,2/7 ,,3/7.-1 



(20) 



= -4.9 x lO-^rads-^^ ^^M^'/^ 1 (21) 

Using Equations (12) and (18) one can find the spin- 
down in wind and field decay solutions 



M (t) 5/7 = m 5 / 7 -5/7C'i 
fi(f) 5/7 = 1 5 / 7 -5/7Ci 



(22) 
(23) 



C 



,_^ mH 7 ns {gm ns ?i 7 
fii 1 



Since the magnetic moment does not change after the 
wind spin-down epoch, Hi at the end of the wind spin- 
down epoch is just the magnetic moment in the present 
RLOF spin-up epoch as estimated above. By assuming 
that f2(<sd) = ^2 << ^1 and Hi « ^1, we obtain 



tsd = - 



7 fti 



5 ^ 2/r M^ s (GM NS yr< 



-3xl0 7 yrh 5 ^% 7 M-t/l 13 M 1A 

(24) 

where we have assumed that /xi = fi as found at the 
end of the previous section. 

From the calculations presented in Equation (24) we 
see that for a captured wind of M W i n d,NS ~ 10~ 2 ■ 
M W i n d ~ 10 13 — 10 14 gs _1 , the spin-down epoch lasts for 
a time 10 7 ~ 8 yr, much shorter than the age of the cluster 
(~ 10 10 yr) for an assumed Hi = Ho = 10 29 ~ 30 Gcm 3 . If 
more wind is accreted, as described for example in the 
focused accretion model of Podsiadlowski & Mohamed 
(2007) , the wind spin-down epoch would be further 
shortened if the wind still carries low angular momen- 
tum. If the wind has instead large angular momentum 
because it is focused through the Lagrangian point L\, 
the RLOF epoch would be further shortened. 

4.4. Incompatibility with a Spin Equilibrium Scenario 

In the previous section we have seen that IGR J17480- 
2446 is not close to equilibrium since it exhibits a strong 
spin- up during the outburst. Several other considera- 
tions on the evolutionary history of IGR J17480-2446 
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also suggest that the NS is currently spinning up. Ac- 
cording to the spin evolution history, the present 11 Hz 
spin frequency does not represent a special frequency. It 
remains to be explained whether the conclusion that the 
binary has currently completed only the first few million 
years since the onset of the RLOF, a phase which lasts 
for ~ 10 9 yr, requires any special pleading. 

If we assume that the observed spin-up and outburst 
duration represent the typical values of an outburst of 
IGR J17480-2446 , then the pulsar spins up by v A t a ~ 
7 x 10~ 6 Hz at each outburst, where A t a is the outburst 
length (55 days). If the NS is dominated by magneto 
dipole spin-down during quiescence, as proposed for some 
accreting millisecond pulsars (Hartman et al. 2008; Hart- 
man et al. 2009; Hartman et al. 2011, Patruno et al. 2010, 
Papitto et al. 2011, Riggio et al. 2011) then a balance of 
the spin-up requires a magnetic field B ~ 10 11 G and a 
recurrence time of ~ 15 yr to provide a spin-down with 
the necessary magnitude calculated via the expression 
(Spitkovsky 2006): 



_ M 2 (¥) 3 (l+smV) 



(25) 



where a is the offset angle between the rotational and 
magnetic field axes. The spin down would then be of the 
order of -10" 14 Hz s" 1 . 

However, as discussed in Section 1, different estimates 
of the B field of IGR J17480-2446 suggest B ~ 10 9 - 
10 10 G. Such a B value requires a recurrence time of 
~ 10 2 — 10 3 yr which is not compatible with a binary 
evolution scenario with the donor in an RLOF phase. 
The mass transfer rate (equal to the outburst mass ac- 
cretion rate averaged over the source duty cycle) would 
then average to M ~ 10 13 — 10 14 gs _1 , a value too small 
for an evolved donor or even a main sequence star with 
M > 0.8 Mq (typical values being - 10 15 - lO^gs" 1 ). 

5. CONSTRAINTS ON THE DONOR MASS AND BINARY 
PARAMETERS 

The mass function of IGR J17480-2446 has been re- 
ported in Papitto et al. (2011) and corresponds to a min- 
imum donor mass of O.4M for a neutron star of 1.4 M . 
The location of IGR J17480-2446 in the globular cluster 
Terzan 5 allows to place further constraints on the donor 
mass. Terzan 5 is composed by two populations of stars: 
one with sub-solar metallicity (Z = 0.01, Y = 0.26) and 
one with supra-solar metallicity [Z = 0.03, Y = 0.29; 
see Ferraro et al. 2009). One interpretation of these dis- 
crepant metallicities and Helium abundances is that the 
cluster is composed by two different populations of stars 
with ages of 12 ±1 Gyr (metal poor) and 6±2 Gyr (metal 
rich). If we assume that the donor star belongs to the 
first population, strong constraints can be placed on the 
donor star. Since the orbital period of IGR J17480-2446 
is « 21 hr, and assuming a minimum NS mass of 1.2 Mq, 
the minimum orbital separation of the system is A w 0.02 
AU. By using the Roche lobe approximation of Eggleton 
(1983): 

0.49 A q^ 
L 0.6<z 2 / 3 + In (1 + qV3) ^> 

with q = Md/M^s being the ratio of the donor and ac- 
cretor mass, the minimum possible Roche lobe radius 



that corresponds to the minimum orbital separation and 
the minimum mass ratio (M d = 0.4 Mq, Mns = 1.2 Mq) 
is l.li? . After 12 Gyr all stars with a mass M>1M 
have evolved to the red giant phase and might have al- 
ready turned into white dwarfs whereas all stars with 
M^0.9M Q (i.e., the turn off mass of the metal-poor 
population of Terzan 5, Ferraro et al. 2009) are still on 
the main sequence. Considering the other extreme case 
of M d = 1.0 M and M NS = 2.0 Mq, the maximum 
Roche lobe radius will be 1.8i? . 

To fill the Roche lobe therefore any possible donor star 
of IGR J 17480-2446 has to be an evolved star that has 
left the main sequence and has increased its radius to fit 
its Roche lobe. The donor star mass falls in the range 
0.9 < M< 1 M and is likely to be a sub-giant. This 
result strongly constrains the orbital separation of the 
system to be in the range 0.023-0.026 AU for a total mass 
of the binary between 2.1 (1.2 M Q NS and 0.9 M Q donor) 
and 3 M (2 M NS and 1 M donor). The inclination i 
of the binary is then constrained by using the projected 
semi-major axis of IGR J17480-2446 (see Table 1): 



sin 



5 x 10~ 3 AU M NS + M d 
A W d 



30° 



(27) 



Identical considerations apply if the donor star belongs 
to the metal rich population, with the difference that the 
turn off mass will be higher by a few tenths of solar mass 
and the inclination smaller by a few degrees. 

All these considerations assume that irradiation of the 
donor star is unimportant in determining its radius. This 
is a good assumption given that non-negligible irradia- 
tion sets in with the onset of RLOF and we have seen 
that this phase has lasted for ^ 10 7 — 10 8 yr. Therefore 
any possible donor star has filled the Roche lobe follow- 
ing a standard evolutionary path that does not involve 
irradiation of the outer envelope. The effect of irradia- 
tion might instead be important during the RLOF phase 
when the outermost layers of the envelope expand after 
absorbing the energy deposited by the X-ray radiation. 
The star being constrained to fit the Roche lobe will in- 
crease its mass loss through the inner Lagrangian point 
and boost the mass transfer rate. The average X-ray 
power absorbed by the companion during an outburst is 



J abs 



(28) 



with (L x ) ~ 4 x 10 37 ergs _1 being the average X-ray 
luminosity during an outburst. Given a duty cycle of 
A 0.005, the long term average power absorbed by the 
donor will be 55 10 33 — 10 34 ergs _1 . This value corre- 
sponds to ^ 10% — 100% of the outgoing power produced 
by a low mass stars in the sub-giant branch and cannot 
be considered negligible. The quantitative effect of irra- 
diation can be calculated with stellar evolutionary codes 
and is beyond the scope of this work. 

6. DISCUSSION 

We have investigated the evolutionary history of 
IGR J 17480-2446 discussing three distinct evolutionary 
epochs: a dipole dominated spin-down epoch, a wind 
epoch and the currently observed RLOF phase. By tak- 
ing into account the effect of timing noise in the X-ray 
pulsations we are able to confirm the long-term average 
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spin-up of 1.4 x 10~ 12 Hzs _1 observed during the 2010 
outburst (Papitto et al. 2011; Cavecchi et al. 2011). If 
the X-ray flux represents most of the bolometric flux, 
the short term spin frequency derivatives v do not scale 
with the X-ray flux in the expected way. The long- 
term z>, however, leads to exclude a scenario in which 
IGR J17480-2446 has already reached spin equilibrium. 
A long-term spin equilibrium with magnetic dipole spin 
down in quiescence balancing the accretion torques in 
outburst is also excluded based on binary evolution con- 
siderations. The spin-up timescale suggests that IGR 
J17480-2446 has been spun-up in the RLOF phase for 
a few 10 7 yr. An interpretation of the observed QPOs 
as the equilibrium frequencies of IGR J 17480-2446 indi- 
cates a similar timescale for the spin-up epoch up to the 
present. 

The spin-down epochs suggest that the NS in IGR 
J17480-2446 had reached a spin frequency of the order 
of ~ 1Hz at the onset of the RLOF phase, if we assume 
that the newborn NSs had typical initial values for fi 
and hq. This reinforces the idea that the time elapsed 
in RLOF has been very short, since the time required 
to spin- up a 1 Hz NS to 11 Hz with a spin frequency 
derivative of 10~ 12 Hzs _1 and a duty cycle of 0.01 is a 
few 10 7 yr. 

Based on these findings we conclude that IGR J17480- 
2446 is in an exceptionally early RLOF phase. The total 
spin- up timescale to transform IGR J 17480 -2446 from a 
slow pulsar (<~ 1 Hz) into a millisecond one (y > 100 Hz) 
is Z 10 8 — 10 9 yr. Since the current RLOF phase will last 
for <~ 10 9 yr, we would expect to observe today ~ 1 — 10 
accreting millisecond pulsars that have followed a similar 
evolutionary history as IGR J 17480-2446 . The accreting 
millisecond pulsar SAX J1748.9-2021 (v ~ 442 Hz) in the 
globular cluster NGC 6440 has a companion with mass 
and orbital period compatible with being in a slightly 
evolved post main sequence phase (Altamirano et al. 
2008). It is likely that the companion of SAX J1748.9- 
2021 has followed an evolutionary history similar to IGR 
J 17480-2446 with the difference that it spent a longer 
time in the RLOF phase thus turning its NS into a mil- 
lisecond pulsar. The scenario outlined above is compat- 
ible with the observed population of accreting pulsars, 
although it is difficult to assess the problem in a robust 
way given the small number statistics. It is also possible 
that dynamical interactions have played a role by de- 
creasing the number of similar systems in globular clus- 
ters via ionization interactions. This however requires 
further investigations along with detailed binary evolu- 
tion calculations. 

The prior spin-down history gives results which are not 
compatible with a binary having an age of several 10 9 yr 
as expected if IGR J 17480-2446 is primordial or if it has 
not suffered exchange interactions in the last few 10 9 yr. 
Different formation scenarios are available to explain the 
apparent discrepancy between the age of the cluster and 
the binary. A first possibility is that a recent dynamical 
encounter has played a role in forming the binary or in 
accelerating the onset of the RLOF phase. If an exchange 



interaction has taken place in the last few 10 7 — 10 8 yr 
then this would explain why the apparent age of IGR 
J 17480-2446 is so short compared with the age of the 
cluster. The location of IGR J17480-2446 in the high 
mass and high central density cluster Terzan 5 suggests 
that the interaction rate might be particularly high in 
this environment. Dynamical simulations to calculate 
the typical interaction rates in the core of Terzan 5 are 
required to assess this question. 

A second possibility for the origin of IGR J17480-2446 
is formation of the accreting pulsar via accretion induced 
collapse (AIC) of a massive white dwarf (Miyaji et al. 
1980, Nomoto 1987). In this scenario a ONeMg white 
dwarf accretes matter until its mass exceeds the Chan- 
drasekhar limit and it collapses to form an NS via elec- 
tron captures on Mg and Ne nuclei. In this case the 
binary must have gone through a preliminary contact 
phase during which the white dwarf was accreting from 
a donor star in RLOF. During the collapse approximately 
0.2 Mq are ejected from the binary (see for example van 
den Heuvel 2011 for a discussion) which causes a sudden 
expansion of the orbit turning the system into a detached 
binary. At this point the binary follows the three evolu- 
tionary epochs described in Section 4 with the NS age set 
by the onset of AIC. The formation of NSs from AIC has 
been discussed by Lyne et al. (1996) to explain the young 
radio pulsar population in globular clusters. Lyne et al. 
note that the formation rate of young pulsars in globu- 
lar clusters is larger than that of millisecond pulsars in 
globular clusters. The young pulsars are found in metal 
rich globular clusters with high mass and high central 
densities, like Terzan 5 (Lyne et al. 1996; Boyles et al. 
2011). Further evidence that globular clusters might con- 
tain young NSs was also presented by Freire et al. (2011), 
who reported the discovery of the youngest millisecond 
pulsar ever found in the globular cluster NGC 6624, with 
a characteristic age t < 24 Myr (see also Lyne et al. 1987 
and Cognrd et al. 1996 on the "young" millisecond pul- 
sar PSR 1821-24A in the globular cluster M28). We can 
therefore speculate that a similar evolutionary process 
has taken place for IGR J17480-2446 thus explaining the 
relatively short lifetime of the NS through the dipole and 
wind spin-down and present RLOF spin-up epochs. 

A possible test for the AIC scenario can be performed 
by measuring the gravitational mass of the NS in IGR 
J 17480-2446 . Since the RLOF phase has lasted for a 
very short time, the NS mass is currently almost identical 
to the value expected from the AIC scenario, which is 
M ~ 1.25 Mq. A gravitational NS mass larger than this 
value would immediately rule out the AIC scenario. 
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